THE QUANTUM DYNAMICS OF THE COMPACTIFIED 
TRIGONOMETRIC RUIJSENAARS-SCHNEIDER MODEL 



J. F. VAN DIEJEN AND L. VINET 

Abstract. We quantize a compactified version of the trigonometric Ruijsenaars- 
Schneider particle model with a phase space that is symplectomorphic to the com- 
plex projective space CP^. The quantum Hamiltonian is realized as a discrete 
difference operator acting in a finite-dimensional Hilbert space of complex func- 
tions with support in a finite uniform lattice over a convex polytope (viz., a re- 
stricted Weyl alcove with walls having a thickness proportional to the coupling 
parameter). We solve the corresponding finite-dimensional (bispectral) eigenvalue 
problem in terms of discretized Macdonald polynomials with q (and t) on the unit 
circle. The normalization of the wave functions is determined using a terminating 
version of a recent summation formula due to Aomoto, Ito and Macdonald. The 
resulting eigenfunction transform determines a discrete Fourier-type involution in 
the Hilbert space of lattice functions. This is in correspondence with Ruijsenaars' 
observation that — at the classical level — the action-angle transformation defines an 
(anti)symplectic involution of CP^. From the perspective of algebraic combina- 
torics, our results give rise to a novel system of bilinear summation identities for 
the Macdonald symmetric functions. 
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1. Introduction 

In a recently published work [R3], Ruijsenaars presented a detailed study of the 
dynamics of the classical Sutherland-Moser particle model [Mo] and its "relativistic" 
deformation the trigonometric Ruijsenaars-Schneider model [RS, Rl]. In addition, he 
also considered a closely related integrable system characterized by an (A^+l)-particle 
Hamiltonian of the form 

(1.1, H= E n ('- 

The Hamiltonian in (1.1) differs from the standard trigonometric Ruijsenaars-Schneider 
Hamiltonian by the substitution (5 ^ ifS (where i — \f—l). Even though over the 
complex field both systems are equivalent, it turns out that their real (i.e. physical) 
dynamics are quite distinct. (Throughout we arc assuming that our variables Xj, 
Pj as well as the scale factors a, (5 and the couphng parameter g are real- valued.) 
The main point is that H (1.1) is periodic not only in the x but also in the p vari- 
ables. This periodicity naturally prompts one to employ a phase space which — upon 
restricting attention to the relative motion in the center-of-mass frame — is bounded 
(in fact compact after a suitable completion). This is in contrast to the situation for 
the standard trigonometric Ruijsenaars-Schneider model (with cos(/3pj) cosh(/3pj) 
and — s\v?{a(3g /2) +sivA^[al3g/2) substituted in (1.1)), where the phase space is 
given by the (manifestly noncompact) cotangent bundle over the configuration space. 

From now on the system determined by the Hamiltonian H (1.1) will be referred 
to as the compactified trigonometric Ruijsenaars-Schneider model. It is the purpose 
of the present paper to investigate the corresponding quantum system. We will see 
that, in accordance with physical intuition, the Hilbert space for the quantum model 
becomes finite-dimensional. In essence, the Hilbert space in question consists of the 
space of all complex functions with support in a finite uniform lattice (grid) over 
classical configuration space. This configuration space has the geometry of a convex 
polytope consisting of a restricted Weyl alcove with walls that have a thickness deter- 
mined by the value of the coupling parameter g. Matching the lattice so as to let it 
fit precisely over the configuration space, including the vertices (corner points) of the 
polytope, produces a quantization condition on g that relates the coupling parameter 
to the size of the lattice. The quantum Hamiltonian is in turn given by a discrete 
difference operator with a step size that is equal to the distance between neighbor- 
ing lattice points. Mathematically, our quantization condition on g translates in 
vanishing conditions for the coefficients at the boundary lattice points, therewith 
guaranteeing that the discrete difference operator Hamiltonian is well-defined and 
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self-adjoint as an operator in the Hilbert space of complex functions over the finite 
lattice. 

For the quantum version of the standard trigonometric Ruijsenaars-Schneider model 
it is well-known (see e.g. the introduction of [D2] or Sect. 7.6.2 of [R4]) that the 
eigenfunctions may be expressed as a product of a factorized (ground-state) wave 
function and Macdonald polynomials (with < g < 1) [M2, M3, M4]. Here also, 
in the case of the compactified Ruijsenaars-Schneider model, the eigenfunctions turn 
out to be similarly expressible in terms of Macdonald polynomials. In contrast to 
the standard situation, however, in the compactified/discrete context of the present 
paper the parameters q and t lie on the unit circle and the diagonalization of the 
model involves only a finite number of Macdonald polynomials (viz., precisely as 
many as the number of lattice points = the dimension of the Hilbert space). The 
symmetry relations for the Macdonald polynomials [Kol, M4, EK] have as conse- 
quence that the discrete kernel for the finite-dimensional eigenfunction transform is 
symmetric. This reflects the fact that we arc actually dealing with a multivariate 
finite-dimensional doubly discrete bispectml problem in the sense of Duistcrmaat and 
Griinbaum [DG, W, G]. More concretely, the discrete eigenfunction kernel satisfies 
the same discrete difference equations in the "spectral" variables as it does in the 
"spatial" variables. Combined with the unitarity, the symmetry of the kernel fur- 
thermore implies that the eigenfunction transform determines a discrete Fourier-type 
involution in the Hilbert space of lattice functions. This is the quantum counterpart 
of the corresponding property of (the closure of) the action-angle transformation for 
the classical compactified Ruijsenaars-Schneider model, which turns out to define an 
involutive (anti)symplectomorphism of the classical phase space (= CP'^) [R3]. 

The paper is organized as follows. We first recall in Sect. 2 some properties of the 
classical compactified Ruijsenaars-Schneider system taken from [R3]. Specifically, we 
discuss the commuting integrals, the configuration space (viz. the restricted Weyl 
alcove with walls of thickness proportional to g) and also the phase space of the model. 
A rather remarkable property of the dynamical system under consideration is that 
the phase space for the relative particle motion in the ccntcr-of-mass frame becomes, 
after a suitable compactification, isomorphic to the complex projective space CP'^ . In 
particular, globally the compactified phase space does not have a topology of product 
form (it is not topologically equivalent to the direct product of the configuration space 
times a (real) A^-dimensional torus). Sect. 3 goes on to demonstrate how canonical 
quantization {pj — > d/idxj) gives rise to discrete difference operator Hamiltonians 
acting in a finite-dimensional Hilbert space of lattice functions over the classical 
configuration space. In Sect. 4, the spectrum of the quantum model is determined 
in explicit form and the corresponding wave functions are expressed in terms of 
Macdonald polynomials with q, t on the unit circle. In order to normalize our wave 
functions such that their L^-norms are equal to one, it is necessary to evaluate a 
terminating version of a recently found summation formula due to Aomoto, Ito and 
Macdonald [Ao, I, M5]. The details explaining how to truncate the Aomoto- Ito- 
Macdonald sum so as to arrive at its terminating version are relegated to the first 
of two appendices at the end of the paper (Appendix A). In a second appendix 
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(viz. Appendix B), some useful properties of the Macdonald symmetric functions 
taken from [M2, M4] have been collected. These properties were needed in Sect. 4 for 
the diagonalization of the quantum model. We have also taken the opportunity to 
reformulate here some of our results from the viewpoint of algebraic combinatorics. 
This leads us, in particular, to a new system of bilinear summation identities for the 
Macdonald symmetric functions (cf. Proposition B.2). The paper closes in Sect. 5 
with some miscellaneous results and remarks. Among other things, it is pointed 
out that: (i) the results on the eigenfunctions give rise to a Discrete Fourier- type 
Transform for lattice functions over the restricted Weyl alcove, (ii) both the maximal 
and the minimal energy of the compactified Ruijsenaars-Schneider model are at the 
quantum level the same as at the classical level (the quantization discretizes the 
energy levels but does not shift the spectrum) and (iii) the dimension of our Hilbert 
space is in agreement with the dimension predicted by the Riemann-Roch-Hirzebruch 
formula for CP^, in the framework of geometric quantization [HK, Si, Hu]. 

Note: By means of a symplcctic (with respect to the standard symplcctic form 
^jdxj A dpj) rescaling (x, p) (/3x, /3~^p) one absorbs the scale parameter j3 in a 
(cf. (l.l)). In this paper we will from now on pick (3=1 without loss of generality. 
At the quantum level this means that we have scaled our variables such that the 
step size of the discrete difference operator Hamiltonians becomes equal to one (or 

to ^ after projection onto the center-of-mass hyperplane), cf. Sect. 3.1 and the 
remark at the end of Sect. 3. 

2. The classical system 

This section serves to summarize some of the basic properties of the classical com- 
pactified Ruijsenaars-Schneider model that are discussed in more detail in [R3]. Since 
we are primarily interested in the relative particle motion in the center-of-mass frame, 
it will be convenient to employ root system notation. For our purposes it suffices to 

restrict attention to the root system of type An- Some relevant preliminaries have 
been collected in the first subsection. For further information regarding root systems 
the reader is referred e.g. to (the Planches in) Bourbaki [B]. 

2.1. Some notational preliminaries. Below the vectors ei, . . . ,ejv+i always rep- 
resent the unit vectors constituting the standard basis of ]R^+-^ and (•, •) denotes the 
(usual) inner product with respect to which this standard basis becomes an orthonor- 
mal basis (i.e. (6^,8^) = Sj^k)- 

Let E be the center-of-mass hyperplane 

(2.1) £; = {xeM^+^ |xi + --- + xjv+i = o}. 

A natural basis {ai, . . . , a.]y} for E is given by the simple roots 



(2.2) 



= - Gj+i, J = 1, . . . ,N. 
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The associated dual basis {cui, . . . ,00^} — determined uniquely by the property that 
{cuj, SLk) — 5j^k — is realized explicitly by the fundamental weights 

(2.3) ujj = {e, + --- + ej) - ;^^(ei + • • • + e^v+i), j = l,...,N. 

To these two bases of (2.1) one can associate the root lattice 

(2.4) Q = Z-Span{ai,... ,aiv} 
and the weight lattice 

(2.5) A = Z — Span{a;i, . . . , uj^} 

as well as the corresponding positive semi-lattices (or integral cones) 

(2.6) g+ = N-Span{ai,... ,aiv} 
and 

(2.7) A+ = N - Span{a;i, . . . , un}, 

respectively (where in our conventions the set of natural numbers N does include the 
number zero). The semi-lattice A+ (2.7) is usually referred to as the cone of dominant 
weights. This cone is partially ordered by the dominance order, which is defined for 
A, /X e A+ by 



+ 



(2.8) n^X iff X-fieQ 

(and 11 -< X iS 11 ^ X and X). 

The Weyl group generated by the reflections in planes orthogonal to the simple 
roots ai, . . . , ajv (2.2) is realized explicitly as the group of permutations a e Sn+i 
acting on the vectors ei, . . . , Cat+i by 

(2.9) (7(ej) :=e^(j). 

The (unique) orbit of the basis vectors a^ with respect to the ^'jv+i-action consists of 
the roots 

(2.10) Ai, = {ej-ek\l<3^k<N + 1). 
For future reference we also need to identify the positive roots 

(2.11) Al^ = An r\Q+ = [Bj - ek\l < j < k < N + 1}, 
the maximal root 

(2.12) ?L^ax = ^ a^- = ei - e^+i ^Ui^un 

l<j<N 
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(this root is maximal in (2-10) with respect to the partial dominance order in 
(2.8)), and the weighted half sum over the positive roots 

(2.13) ^ - f E ^ = f E (^^- - 

aeA+ l<j<k<N+l 

= I ^ (iV-2(j-l))e, 

l<j<N+l 

= g{uji^ \-ujn). 

2.2. Integrability. The Hamiltonian H (1.1) is known to be integrable: a complete 
set of integrals in involution is given explicitly by [RS, Rl] 

(2.14) j: cos(E,,,p,)n(i- ^,^.„ /.l' . 

\J\=r 

r — 1,... ,A^ + 1. (Recall that we have rescaled the variables such that (3 — 1, cl. 

the note at the end of the introduction.) Observe that Hr (2.14) specializes for r = 1 
to the Hamiltonian H (1.1) and that H^^i = cos(pi + ■■■ +pN^i), reflecting the 
translational invariance of the model. The projection of the Hr-Aow onto the center- 
of-mass hyperplane xi + ■ — h xn+i = is governed by the reduced Hamiltonian Hr, 
written conveniently in root system notation as 



(2.15) Hr^ cos((i/,p)) H 



(a,i.)=l 



sin^(f) 



sin f(a,x) 



r — 1, . . . ,N (for i^AT+i the reduced flow in center-of-mass plane is of course trivially 
stationary). Here x := {xi, . . . ,xi\f^i), p := {pi, . . . ,pAr+i) and the sum in (2.15) is 
over all weights i/ e A (2.5) that lie in the 5'Ar+i-orbit (recall the action (2.9)) of the 
r-th fundamental weight vector cj,. (2.3). 

2.3. The reduced phase space for the relative particle motion: CP^. Let us 
from now on assume that the scale factor a is positive and that the parameter g lies 
in the interval 

In order to arrive at real- valued Hamiltonians Tir (2.15), one is led to employ a 
conflguration space in which the particle distances \xj — Xk\ are bounded from below 
by gf (> 0) and from above by 27r/a — g (> 0). This is realized by picking as 
conflguration space the submanifold Tig of the center-of-mass plane consisting of the 
points X G (2.1) satisfying the conditions 

(i) (aj,x) > ^ for j = 1,... ,Ar; 

(ii) (a^aa;,x) < 271 / a - g 

(where the vectors ai, . . . ,aN denote the simple roots (2.2) and Simax is the maximal 
root (2.12)). The parameter restriction (2.16) ensures that the submanifold T,g C E 
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determined by (i), (ii) is nonempty (add the inequalities from (i) and use (2.12) to 

compare with (ii)). Furthermore, has the geometry of an open convex polytope 
consisting of an alcove with walls of thickness inside the Weyl alcove Sq (which 

corresponds to the limit g [ Q). The open convex polytope (or open simplex) is 
completely determined by the A'^+l vertices (corner points) p, p+Mujr (r = 1, . . . ,N) 
with M = ^ - (AT + 1)^ > 0. See Figure 1. 



Figure 1 . The restricted Weyl alcove with walls of thickness g/ -\/2 for 
N — 2. The region of the inner alcove corresponds to the configurations 
spaces (without boundary) and Hg (with boundary) of the three- 
particle system in the center-of-mass hyperplane Xi + X2 + = 0. 
The vertices and boundary segments with/without the shifts between 
square brackets refer to the inner/outer triangle, respectively. 

An obvious candidate for the phase space would now of course be the cotangent 
bundle over the configuration space: T*(T,g) ^ T,g x E. However, in view of the 
periodicity of the the Hamiltonians Tir (2.15) with respect to translations in p over 
vectors in the dilated root lattice 27rQ (cf. (2.4)), it is natural to restrict to a smaller 
phase space of the form x T, where T is the A^-dimensional torus E/{27tQ). This 
torus can be coordinatized explicitly as 

(2.17) T={peE\-n< {ur, p) < tt, r = 1, . . . , A^}, 

where the components (0;^, p) of the vector p with respect to the basis of fundamental 
weights {cui, ... ,ujn} should be read modulo 27r. 

Unfortunately, it turns out that the Ti^-flows are not complete on the bounded 
phase space x T [R3]. To remedy this incompleteness, it is needed to compactify 
the phase space in a suitable manner. For this purpose the key observation from 
[R3] is that it turns out possible to embed the noncomplete phase space x T 
densely and symplectically in CP^. Here the complex projective space is to be 
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viewed as a 2A'"-dimensional real manifold with symplectic form proportional to the 
standard symplectic form inherited from the Fubini-Study Kahlcr metric on CP^^. 
The Hamiltonians Tii, . . . ,T-Cn (2.15) lift under this embedding to smooth (Poisson 
commuting) Hamiltonians on CP^ [R3] and the completeness of the corresponding 
Hamiltonian flows is thus immediate from the compactness of the extended phase 
space CP^. 

The relevant embedding of the noncomplete phase space x T into CP^ presented 
by Ruijsenaars is given explicitly by (x, p) i— >• [1 : : 2:2 : • • • : Zn] G CP^ with 

The inverse mapping [zq : Zi : Z2 : ■ ■ ■ : z^] ^ (x, p) G x T, defined on an open 
dense patch {[zq : ■ ■ ■ : zn] \ ^ {j = 0, . . . , iV) } of CP^, reads 







2 




|2 + 




' + 


\zn\ 


2 



(2.19a) (a,-,x) = (27r/a - (TV + 1)^) , ^ 12 , , u 12 + 9, 

(2.19b) e'^^^'P^ -- 



j — 1,. . . ,N. (This gives the components of x, p with respect to the bases {ai, . . . , an} 
and {cui, . . . , con}, respectively.) The above mappings are symplectic when x T is 
equipped with the standard symplectic form induced by dxj A dpj and CP^ is 

endowed with the renormalized Fubini-Study symplectic form 

(2.20) ujR = J J2 ^ - ^iv I 12 {Yl ^^^^i) ^ (Zl ^i^'^i 

l^j=o\^3\ \j=l Z^j=0Pjl j=l j=l 

where the normalization is such that the integral of cur over a complex projective line 
equals AnR"^ with 

(2.21) 2R^ = — -(N + l)g. 

a 

The coordinate functions in (2.19a) clearly extend to smooth functions on the 
whole of CP^. The image of the extension of the coordinate map to the completed 
phase space CP^ is therefore given by the compactification Hg of Eg in E 

— 271 

(2.22) Eg = {x e E I (a,-, x) > (? (j = 1, . . . , iV); (a„„„ x) < — - ^ }. 

It is natural to interpret the simplex E^ as the configuration space for the compacti- 
fied trigonometric Ruijsenaars-Schneider model, even though globally the completed 
phase space CP^ has not a topology of product form. (In particular CP^ ^ Ep x T.) 
Notice in this connection that the coordinate functions in (2.19b) for the momentum- 
like variables do not extend continuously to the boundary hyperplanes zj = 0, 
j = 0, . . . ,N (as the limiting value of (2.19b) for Zj along a radius in the 
complex plane depends on the argument). 

It is quite instructive to view how the compactification works topologically in the 
situation of two particles {N — 1). In this special case the reduced phase space T,gXT 
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before completion has the structure of an open hne segment (Eg = ~^/2) | x e 

2Ti/a — g[ }) times a real one-dimensional torus T^. Topologically this is a cylinder 
without the two boundary circles or, equivalently, a two-sphere with two distinct 
points extracted. The compactification adds the two extracted points (pinching) thus 
resulting in a compact phase space with the topology of a two-sphere S"^ (^ CP^). 
The canonical projection 11 : x T i— > Eg clearly extends uniquely to a continuous 
projection H of S"^ onto the closed line segment E^ = {(2^/2, —x/2) \ x & [g ,271 / a — 
g\ }), however, the fiber IT (m) with m G Eg reduces to a point when m lies on the 
boundary Eg \ Eg whereas it is isomorphic to a real one-torus for m in the interior 
Eg. In particular, we do not have that S"^ is isomorphic to Eg x T^. See Figure 2. 



Figure 2. The compactification of the phase space for = 1 turning 
the cyhnder Eg xT into a two-sphere (=CP^) by adding two points. The 
arrows indicate the canonical projections 11 , 11 onto the configurations 
spaces Eg (open line segment) and Eg (closed line segment), respec- 
tively. 



3. Quantization 

In this section we quantize the compactified Ruijsenaars-Schneider model of the 
previous section. The quantum versions of the Hamiltonians . . . ,1-iN (2.15) for 
the relative particle motion in the center-of-mass frame will be given by commuting 
discrete difference operators acting in a finite-dimensional Hilbert space of functions 
with support on a finite uniform lattice over the classical compactified configuration 
space Eg (2.22). 

3.1. Ruijsenaars difference operators. In [Rl] Ruijsenaars showed how the Pois- 
son-commuting Hamiltonians Hr (2.14) may be quantized formally (i.e., without 
specifying a Hilbert space) by means of canonical quantization in such a way that 
integrability is preserved. For the reduced integrals T-Lj. (2.15) the procedure leads to 
difference operators of the form 

(3.1) nr^{n^ + n;)/2, r = i,...,7v, 
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with 

(3.2a) nt = Yl K'/'(x) V;'/'(-x), 

(3.2b) n; = Yl yy'i-^)T-'vy'{^), 

where 

(3.3) = e<"'^> 

denotes the operator acting on functions f : E ^ Chj a translation over the weight 
vector u, i.e. (Ti^/)(x) = /(x + i/), and the coefficients are determined by 

(3.4) K(x)= n 

sinf(a, x) 

(a,^>=l 

The commutativity of the above difference operators is by no means evident from 
their exphcit expressions, but it does follow immediately from Ruijsenaars' results 
in [Rl]. To this end it is helpful to observe that — 'H%j^i_^ because —ujr is in 
the S'AT+i-orbit of cu^v+i-r and Vt,(— x) = K_j^(x). Hence, the commutativity already 
follows from the commutativity of TiJ*", . . . , 7i^, which can be traced back to [Rl] by 
noticing that H+ corresponds in Ruijsenaars' notation to the operator S^. S^'_l.[^^^^ 
(with the number of particles of course being equal to iV + 1). 

Proposition 3.1 (Quantum integrability [Rl]). The difference operators H^, 7Y~ 
and Tir, r — 1, . . . ,N mutually commute. 

To see that the classical version of Ti.,. indeed amounts to Tir (2.15), one observes 
that after substituting T„ = exp(i(z/, p)) (which is the classical analog of (3.3)) in Tir 
one arrives at 7ir (2.15) by using the identity 

y„(x)y„(-x)= n 

\ sm i7(a, x)/ 



3.2. The finite-dimensional Hilbert space. The formal difference operators Hr 

(3.1) and "Hf (3.2a), (3.2b) shift function arguments over vectors in the weight lat- 
tice A (2.5). We will now assign a precise meaning to these difference operators as 
operators in a Hilbert space of functions with support in a uniform lattice over the 
classical compactified configuration space T,g (2.22). 

The point p (2.13) denotes the "minimal" vertex of the simplex determined 
(uniquely) by the property that the functionals (a^, ■), j = 1, . . . , simultaneously 
assume their minimum value g. By shifting from p over vectors in the weight lattice 
A (2.5), one generates a uniform lattice in Tig consisting of the points p + p, p & A"*" 
(2.7) with {aLjnax:P + 1^) ^ Ng + (a^a^,,//) < ^ - g. When the (positive) coupling 
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constant g and scale factor a are related by 

(3.5) (iV + l)^ = MeN\{0}, 

a 

then the maximum value — g oi the functional (amax, ■) is assumed on the lattice. 
More to the point, it means that in this situation apart from the "minimal vertex" 
p also the other vertices of the simplex Eg (2.22) (viz. the "maximal vertices" 
p + Mur, r — 1, . . . ,N) lie on the lattice and, hence, that the lattice fits precisely 
over the classical configuration space including its boundary. See Figure 3. 



Figure 3. The lattice p + A^^^ over the configuration space sup- 
porting the quantum wave functions for = 2 and Af = 5. The 
dimension of the Hilbert space L^{p + A^) (i.e. the number of points 
in the lattice) amounts in this case to (^2^) =21. 

From now on we will assume that the condition in (3.5) is satisfied. (Notice that 
the condition is compatible with the parameter restriction in (2.16).) Let A^ be the 
alcove of dominant weights in A+ (2.7) given by 

(3.6) A+ = {A e A+ I (a„„„ A) < M] 

and let L^(p + A^y) be the finite-dimensional Hilbert space of complex functions over 
the lattice p -\- h.]^ := { p -\- p \ p & A|^ } endowed with the (standard) sesquilinear 
inner product 

(3.7) {f,h)= f{p + l^)W+l^ 

(/, he L^{p + Ai^)). Notice that our conventions are such that the inner product is 
linear in the first slot and antilinear in the second slot. The dimension of the Hilbert 
space L'^{p + A]^) amounts to the number of points in A^j, i.e., the number of vectors 
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of the form riicui + • • • + tincun with rij & N {j — 1, . . . , N) and {a^ax, nicui + ■ — h 
unoJn) —ni-\ Vun < M: 

(3.8) ata(L^(p.Aa)^M. 

In order to see that the difference operators 'H^'^ (3-1), (3.2a), (3.2b) are weU- 
defined as operators in L'^{p + A^) we need the following lemma. 

Lemma 3.2 (Regularity, positivity and vanishing boundary conditions). For positive 
parameters a,g subject to the condition (3.5), /i G (3.6) and v in the SN+i-orbit 
of a fundamental weight vector ujr (2.3), one has that 

<V^{p + li) < oo, <V^{-p- fj.- u) < oo if /i + z/GA+, 

V,{p + p)=0 if p + u^Klj 
and, if in addition g ^ {-^, jfzi-i ' ' ' ^h^}' ^^^^ 

—oo < Vy{—p — p — u) < oo if p + u ^ A^ 
(where p and K/(x) are given by (2.13) and (3.4), respectively). 
Proof. Let us write 

TT sinf((a,p + /i)+ff) -pp sinf((a,p + /i) -^) 



sin|(a, p + yLt) sin|(a, p + 

(a,i/)=l {a,u)=-l 



Prom the inequality 



^ / X 27r 27r 

(*) 0<g<{sL,p + p)< g< — 

a a 

for a e A'^ and p G A^, it is seen that all factors in the denominator of the 
above formula for T4(p + p) are positive and that all factors in the numerator are 
nonnegative. Zeros in the numerator appear when {a, p + p) + g becomes equal to 
271 /a or when {a, p + p) — g becomes equal to 0. The first situation occurs if and 
only if {amax,^) = 1 and (a„,a^,//) = M, i.e., iff {amax,IJ' + > M . The second 
situation occurs if and only if (a^ , u) = —1 and (a^, p) — for certain simple root a^ 
(2.2), i.e., iff (a^-, p + ly) < for certain simple root a^. 
In a similar way one derives from the formula 

K(-p - p-v) = 

T-T sinf ((a,p + /i) + 1- g) yr sinf ((a,p + ^) - I + g) 
-L-L sinf((a,p + Ai) + l) ^\ sinf((a,p + A^)-l) 

(a,j^) = l (a,i/)=-l 

combined with the inequality (*), that Vv{—p — p — v) is positive and finite for 
p e A^ with p + v & The denominators become zero when (a, p + p) + 

1 — 2T:/a or when (a, p + //) — 1 = 0, which can happen only ii p + v ^ A^ and 

y ^ liV' AT-l' iV-2' ' 2' /■ '-' 
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We learn from Lemma 3.2 that for parameters subject to (3.5) the coefficient 

functions K,(x) (3.4) are regular and positive on the lattice points p + yU, /i G with 
/X + z/ G and zero on the boundary lattice points p + ji, ji E A^ with ji + u ^ K\j. 
Similarly, the coefficient function Vi,(— x — v) is regular and positive on the lattice 
points p+p with p, p+u G A^ and generically (i.e. for g ^ j = 1, . . . ,N) regular 
on the boundary lattice points p + p with p G A^, p + i/ ^ A^. (Shifts of the type 
— > ^x— 1/ in the functions governing the coefficients of (3.2a), (3.2b) originate 
from commuting the coefficients on the right of the translation operator T^^ to the 
left.) This entails that for parameters subject to (3.5) (and g ^ {j^, j^Zi^ • • • ? 1}) the 

difference operators Hi, . . . (3.2a), (3.2b) (and hence Hi, . . . ,Hn (3.1)) admit 
a well-defined restriction to functions in L^(p + Aj^) which maps the Hilbert space 
into itself. The exphcit action on functions / G I/^(p + A^) is given by 

(3.9) {Hff){p + p)^ J2 W^{p + p)f{p + p + u) 

i^eSjv+i('^T-) 

{p G A+ ) with 

u) >0 for p + u e 

for p + u ^ A+ , 

u)>0 for p-veKli 
for p-v ^Klj 



W^ip + p) = \^Vu'\p + p)V^'\-p-p- 
W;{p + p) - 



V^'\-p-p)V:'\p + p 




.1/2, 



(recall to this end also that K(— x) = V-y{yi)). The vanishing boundary conditions 
for the coefficients W^{p + p) guarantee that {Hff){p + p), p e A^ depends only 
on the values of /(•) in the points of the lattice p + A^, i.e., we have that Hf is 
well-defined as operator in L'^{p + A^^). For g E {j^, j^zi^ ■ ■ ■ ) 1} ambiguities in the 
value of the products Vu{p + p)Vu{—p — p — u) and K/(— p — p)Vy{p + p — u) at the 
boundary points p G A^ with p + v ^ A^ may arise as zeros deriving from Vy{p + p) 
and T4(— p — p) can meet with possible singularities deriving from V^(— p — p — u) 
and Vy{p + p — v), respectively. For instance, when = 1 we have for generic x G 
that V^(x)Vt,(— X — v) = 1. Here we will resolve such ambiguities in the value of 
coefficients at the boundary lattice points for g in the exceptional set {jj, ... ,1} 
by requiring continuity under small variations in g. Specifically, this means that for 
all parameter values subject to the condition (3.5) wc will pick the action of H^ with 
vanishing boundary conditions in accordance with (3.9). 

Proposition 3.3 (Self-adjointness). Let us assume (positive) parameters subject to 
the condition (3.5) and let the action ofHf : L'^{p + A'lj) — >• L^(p -|-A^) be given by 
(3.9). Then the difference operatorsHj. — {H^+H~)/2, r — 1, . . . ,N are self-adjoint 
in L2(p + A+). 
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Proof. It suffices to demonstrate that the operators and are each others' 
adjoints in L'^{p + A^). Some elementary manipulations produce: 

E (^rV)(p + /^)%+7^ 

E E /(P + /i)W^."(P + /i)MP + /i-'^) 

i^eSjV+l(wT-) ji6A+ 

/i-i^eA+ 

^'^^■'^ E /(-^ + iiW-h){p + fi) = (/, 7^,-/1), 

AeA+ 

where we have used (i) the vanishing boundary conditions for the coefficients W^, 
(ii) the substitution p — jx — v and (mz) that W^{p + p — u) — W~{p + p) >0. □ 

Remark: For given trigonometric period ^ > 1, the parameter restriction in (3.5) 

determines a quantization condition on the coupling parameter g (permitting only 
a finite number of values for g labeled by M G {1, • • • , [77]})- However, it is also 
possible (and probably somewhat more natural) to instead interpret the restriction 
in (3.5) as a quantization condition on a step size parameter. Recall to this end that 
in the present paper we have scaled our variables such that the scale parameter f3 
appearing in the classical (compactificd) Ruijsenaars-Schneider Hamiltonian H (1.1) 
has the value 1 (sec the note at the end of the introduction). By substituting xj —>■ 
l3~^Xj and a ^ a/3 {/3 > 0) in the difference operators of Sect. 3.1, we reintroduce 
the scale parameter P in our quantum Hamiltonians. Specifically, the operators TYr^^ 
(3.1), (3.2a), (3.2b) then pass over to discrete difference operators of the form given 
in (3.1)-(3.4) with the coupling parameter g and the translation operator (3.3) 
replaced by (3g and T^, = exp(/3(i/, ^)), respectively. In other words, at the quantum 
level the scale parameter P enters as the step size parameter of the discrete difference 
operators [Rl]. As a consequence of this rescahng, the lattice supporting the wave 
functions is going to be scaled by f3 resulting in a lattice of the form /3(p+A^^) and the 
parameter restriction in (3.5) passes over to the condition ^-{N+l)g = M e N\{0}. 

For a given trigonometric period ^ and (positive) coupling parameter g(, the latter 
parameter restriction may be interpreted as a quantization condition on the step size 
parameter /3 (permitting an infinite series of values for /3 labeled by M e N\{0}). The 



(3^9) 
i 
(n) 

{Hi) 
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quantization condition at issue adjusts the step size such that the lattice (5{p + A^) 
fits precisely over the classical configuration space Ylpg (cf. (2.22)) including the 
corner points (vertices). 



4. Wave functions 

In this section an orthonormal basis for L^(p + A^) is presented consisting of joint 
eigenfunctions of TYi, . . . , TY^v- 



4.1. A factorized eigenfunction. Let 

1 



(4-1) A(/.) = ^ , , e A+ 



with 



(4.2a) CM = n 



((a,p) : sin«; 



(^ + : sin„)(a,M) 

^ ^ fi; (l + (a,p):sin,,)(a,^) 

Here we have introduced "trigonometric Pochhammer symbols" defined by 

m = 



(4.3) {z : sin^). 



nr=o^sinf(^ + A;) m = 1,2,3,... 



We need two preparatory lemmas. The first states that, for positive parameters 
subject to (3.5), the value of C±{fi) (and hence that of A(/i)) is positive and finite 
for jj, e A^ (3.6); the second lemma describes a functional relation between A (4.1) 
and the coefficient functions K(x) (3.4), z/ G S]\f+i{ujr)- 

Lemma 4.1 (Regularity and positivity). For positive parameters a, g subject to the 
condition (3.5), one has that 

C±{n) > for A* e A+ 

(with A'Ij given by (3.6)j. 

Proof. Using inequality (*) in the proof of Lemma 3.2, it is not difficult to infer that 
the arguments of the sine factors in C±(/x) (4.2a), (4.2b) lie between and 7i/2 . □ 

Lemma 4.2. Let v be in the SN+i-orbit of a fundamental weight vector ujr (2.3) and 
let 11,11 + V e A+ (2.7). Then 

A(ii + i/)K(-p -ii-u)^ A(//)K(p + A*). 
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Proof. We have that 



(1 + (a, p) : siiiQ,) 
7+ (1 - f + (a, P) : sinc,)(a,^+^> ((a, p) : sin«)(a,^+j.) 

(S' + (a, p) : sina)(a,;.) (1 + (a, p) : sing,) (a,;,) 
^ + (a, p) : sin„)(a,^> ((a, p) : sina)(a,^) 



aeA+ 

X J"!- smf((a,p + p) +^) smf((a,p + /i) + 1) 



X 



sm f ((a, p + p) + smf(a,p + p) 

(a,i/)=l 

sin |((a, p + /i) - ^f) sin|((a, p + p) - 1) 



n 



sinf((a,p + p) - 1 + ^) sin f (a, p + p) 

(a,i/)=-l 



Multiphcation by 



TT sinf((a,p + p) + l-f/) -j-j- sin f ((a, p + p) - 1 + ff) 
A A sinf((a,p + p) + l) ii sinf((a,p + p)-l) 

aeA+ aeyi+ 
{a,i/)=l {a,i/)=-l 

leads after cancellation of common terms in the numerator and the denominator to 
an expression of the form 

{g + (a, p) : sinc,)(a,M> (1 + (a, p) : sin„)(a,^) 



na 



(1 - ^ + (a, p) : sin„)(a,M> ((a, p) : sin«) 



(a,M> 



JV 



^ TT sinf ((a, p + p) +g) -|-|- sin f ((a, p + p) - g) 
sin |(a, p + p) -L-L sinf (a, p + p) 

(a,!/)=l {a,i/)=-l 

= A(p)K(p + /^). 

□ 

After these preliminaries we are now in the position to introduce a factorized joint 
eigenfunction of Tii, . . . , TIn- Let ^ (p + ^ M be the lattice function defined 

by 

(4.4) *o(p + = ^'^'(z^)' e A+ , 

where the normalization constant Mq is chosen such that {^o, '^o) = 1 (recall the 
inner product (3.7)). Notice that is well-defined and positive at the lattice points 
p + p, p e (3.6) because of Lemma 4.1. 
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Proposition 4.3 (Factorized eigenfunction) . For positive parameters a, g subject 
to condition (3.5), the function \l/o (4.4) is a joint eigenfunction of the difference 
operators Hr : L'^{p + A^^) ^L'^{p + Aj^) (3.1), (3.9) 

= ( cosq;(i/,p)) *o, r = l,...,A^ 

(where p is given by (2.13)). 
Proof. One has that 

(7^.*o)(p + /^)^'"^^'-'^-^ E (w^.nP + /^)AV2(p + /^ + -) 



(i) 



(ii) 



+ W-{p + p)A'/\p + p-u)) 

E v.{p + p) + v.{-p-p))a'/\p + p) 

E cosa(i/,p)) ^'o(p + /x), 

i^G5jv+l(t<^r) 



where we have used: (i) Lemma 4.2 combined with the fact that W^{p + p) = 
V±u{p + //) = for // ± 1/ ^ Aj^ and {O) the Macdonald identity [Ml, Theorem (2.8)] 
(see also Appendix B) 

E = E cosa{u,p). 

i^eSjV+l(t^T-) i^6Sjv+l(Wr) 

□ 

The following proposition gives a compact product formula for the proportionality 
constant Afo, which normalizes the wave function (4.4) such that its L^-norm is equal 
to 1. 

Proposition 4.4 (Normahzation). The value of Mo is given by 

Mo=Yl ^(/^) = 2^(^-^)(7V + 1) H {1 + ng: sin,)M-i, 

where it is assumed that the parameters satisfy condition (3.5). 

Proof Clearly A/'o = Z)^eA+ A(/^) normalizes (4.4) such that (^o, ^o) = 1- The 
evaluation of the sum leading to the product formula on the r.h.s. hinges on a 
terminating version of a recent summation formula due to Aomoto, Ito, and Mac- 
donald [Ao, I, M5]. The details of the summation are relegated to Appendix A (see 
(A.4)). □ 

By speciahzing the Macdonald identity in the last line of the proof of Proposition 4.3 
to X = p and recalling Lemma 3.2, one arrives at a simple product formula for the 
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eigenvalues: 

(4.5) Yl ^««"(^'^) = E V.iP)'^~'"VM 

^ nf=?sin(^) 

n;=iSin(^)nf=V"''sin(^)- 
For r = 1 this product formula specializes to the well-known geometric progression 

(4-6) 2^ cos(ap,) = . 

4.2. The complete eigenbasis. We will now extend the factorized wave function 
^0 (4.4) to an orthonormal basis of L^(p + Aj^) consisting of joint eigenfunctions 
of the commuting operators Tii, . . . jUn- The eigenbasis will be expressed in terms 
of Macdonald polynomials with 1^1 = 1. To describe these we need notation for the 

elementary symmetric functions 

(4.7) E^{^) = J2 e^*"<"'"\ r = 1, . . . , AT 
together with their real parts 

(4.8) £^r(x) = cosq;(i/,x), r^l,...,N 

and also for the monomial symmetric functions 

(4.9) mA(x)= J2 e'"^'''"^ AeA+. 

/ue5jv+i(A) 



Notice that E;(x.) = E+{x) = 

i-r(x) (since —ujr lics in the S'Tv+i-orbit of 
ojN+i-r)- The Macdonald polynomials px{x.), A e A"*" are now defined as the unique 
trigonometric polynomials of the form 

(4.10a) pa(x) = mx{x) + ^ cx,u.mf,{x), cx,^, e M 

i[ieA+ 

satisfying the difference equations 

(4.10b) K(x)pa(x + z/) = ^+(p + A)pa(x), r=l,...,Ar 

(with the coefficients K/(x) given by (3.4)). For generic parameters the existence of 
polynomials px of this form follows from the work of Macdonald in [M2, M3, M4] 
(see Appendix B for a brief summary of the results most relevant to us here). Our 
parameters a and g are related to the parameters q and t employed by Macdonald 
via t — and q — e*" (cf. Appendix B). So, in particular, we have that |?| = 1 (and 
\t\ — 1) for real a (and g). 
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An important property of the Macdonald polynomials is that after renormalizing 
in the following way 

(4.11) Pa(x) = C+(A)pa(x), AeA+ 

(where C+(A) is given by (4.2a)), they satisfy the symmetry relations [Kol, M4, EK] 
(cf. also Appendix B) 

(4.12) Py^{p + ^i) = P^{p + \), A,/xeA+. 

The polynomials in (4.11) are normalized such that Px{p) = 1, as is clear from the 
symmetry relation (4.12) specialized to = (since Po{-) = 1). 

Even though for generic parameters the existence of the Macdonald polynomials 
of the form (4.10a), (4.10b) is guaranteed by Macdonald's work, it is a priori not 
entirely obvious that it is possible to specialize them to positive parameter values 
for a,g subject to the constraint in (3.5). The point is that for certain special 
values of the parameters the eigenvalues E^{p-\- A) on the r.h.s. of (4.10b) may not 
be semisimple. This manifests itself through possible singularities in the expansion 
coefficients c^,^ of (4.10a) at such special parameter values. The next lemma ensures 
that for A G A^ the eigenvalue E^[p+\) is in fact semisimple for positive parameters 
a, g subject to condition (3.5) and, hence, that the Macdonald polynomials A e 
Aj^ indeed admit a well-defined specialization to these parameters values (without 
any singularities in the expansion coefficients being hit). 

Lemma 4.5 (Semisimple spectrum). Let a, g he positive and subject to the condition 
(3.5). Then the elementary symmetric functions E^{yL), . . . , E^{'x) (4.8) separate the 
points of the lattice p + A~j^. 

Proof For x,y E E (2.1), we have that -B+(x) = £;+(y) for r = 1, . . . , if and only 
if 

x = a(y) mod —Q, 

a 

with a e Sn+i- If both x and y lie in the Weyl alcove Sq (i.e. the open simplex 
characterized by the conditions (i) and (ii) of Sect 2.3 with g 10), then the only way 
in which this is possible is if cr = id and x = y. (The Weyl alcove Sq determines 
a fundamental domain for the action of the affine Weyl group Sn+i ^ i^Q) on E.) 
The lemma then follows because the conditions on the parameters guarantee that 
p + A+ C Eg C So. □ 

After these preliminaries let us now introduce the wave function : (p+-^m) ~^ *^ 
given by 

(4.13) *,(p + /.) = -^AV2(A) AV2(/.)P,(p + /.), A, e A+ . 

Notice that for A = this wave function reduces to the factorized wave function 
in (4.4). The next symmetry property is an immediate consequence of the symmetry 
relations (4.12) for the renormahzed Macdonald polynomials -Pa(x) (4.11). 
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Proposition 4.6 (Symmetry). One has that 

^A(p + /i) = ^m(p+A) for A,/iGA+. 

The function "i/x (4-13) turns out to be a joint eigenfunction of the operators 
Til, ■ ■ ■ , Ti-N- 

Proposition 4.7 (Diagonahzation). Let us assume positive parameters a,g subject 
to the constraint (3.5) and let A e Aj^. Then 

nr^x^Er{p + X)^x, r^l,...,N 

(where fir and Er{-) are given by (3.1), (3.9) and (4.8), respectively, and p is taken 
from (2.13) ). 

Proof. One has that 

i^e5jv+i(a)T-) 

i -^AV^(A)AV^(/.) X 

V,{p + p)Px{p + p + iy) 

E+{p + X)^lfx{p + p), 



(n) 



where wc have used (i) Lemma 4.2 combined with the vanishing properties of the 
coefficients Wj~ and at the boundary (cf. Lemma 3.2), and (ii) the defining 
difference equations for the Macdonald polynomials in (4.10b). The proposition now 
follows from the observation that Tir = {Ti.^ + 7i~)/2 and that Er{-) = {E^{-) + 
E-i-))/2 with n; = iiti+i-r and E;{-) = £;++!_,(•)• □ 

It is clear from the proof of the above proposition that the functions (4.13) in fact 
diagonalize the operators (3.9) individually 

(4.14) n^^x^E^{p + X)^x, AeA+ 

(with the parameters satisfying (3.5)). 

It is a priori not obvious that the functions ^a, A G A^^ actually span the Hilbert 
space L^(p+ A^^), since in principle linear dependencies might arise between the Mac- 
donald polynomials p\{-), A G A^ upon restriction to the lattice p + A^. However, 
the following result states that the functions ^a, A G A^ in fact form an orthonormal 
basis of L^{p + A^), therewith excluding the possibility of such linear dependencies 
to occur. Phrased alternatively: the lattice evaluation homomorphism from the poly- 
nomial subspace Span{m\} ^^^Jy+ to the space of complex functions over the lattice 

p + A'l^ — defined by the assignment m\{x) i— >• mx{p + p) — is an isomorphism. 
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Proposition 4.8 (Orthonormality). For positive parameters a,g subject to condi- 
tion (3.5), the functions "ifx : (p + A^) — > C m (4.13) form an orthonormal basis of 
L^p + A+j), i.e. 

\ 1 if \ = n 

Proof. The orthogonality follows by applying the eigenvalue equation (4.14) to the 
equality (7iI^\E'a, \E'^) = (\E'a, Ti^^E'^) (cf. the proof of Proposition 3.3) and using that 
the elementary symmetric functions E^{-), . . . , Ejf{-) separate the points of p + 
(cf. Lemma 4.5). 

To see that the normalization of the wave functions is such that their L^-norms 
are equal to 1, we apply the symmetry relations (Proposition 4.6) to the eigenvalue 
equations of Proposition 4.7. This leads to a system of difference equations for the 
wave functions in the spectral variable of the form 

^,(p+A)^^(p + A) = 

E 

r = 1, . . . , (where is taken from (3.9)). Applying the expansion on the r.h.s. 
to both sides of the equahty {Er^^, ^n+ur) — i^n^ -E'r-^/i+wr-) ^iid exploiting the 
orthogonality of the wave functions, produces the relation 

But then, since Wj^(p + p) — W~^{p + p + oUr) (> 0) for p,p + Ur G A^ (see (3.9)), 
it is immediate that (^^, ^^) is independent of /i G A^. The orthonormality now 
follows because ior p — we have that (^o, ^o) = 1 in view of Proposition 4.4. □ 

Remarks: i. The orthonormality of the wave functions (4.13) described by 
Proposition 4.8 can be rewritten in terms of discrete orthogonality relations for the 
Macdonald polynomials Pa(x) (in the monic normalization) or -Pa(x) (in the sym- 
metric normahzation with P\{p) — 1, cf. (4.11), (4.12)). The discrete orthogonality 
measure is supported on the lattice p + A|^ with positive weights given by A (4.1). 
Specifically, we conclude from Proposition 4.7 that for positive parameters a,g sub- 
ject to the condition (3.5) one has that 

tf p^X 



(w^+(p + p)^,+.{p + A) + W-{p + p)^,..{p + A)) , 



(4.15a) Yl P^(P + ^) Pt^(P + ^) ^(^) = i A/- 



or equivalently 



if p = X, 



(4.15b) Yl ^a(p + ^) P,{P + ^) H^) = ( •{ ^.t\ 

L A(A) ''J P — ^' 

for A, p e A^ (with C± and Mq given by (4.2a), (4.2b) and Proposition 4.4, respec- 
tively) . 
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ii. If we associate to each dominant weight vector A given by 

N 

(4.16a) X^^ljUj, IjeN 

a contragredient dominant weight A* of the form 

N 

(4.16b) X* = J21n+i-jOJj, 

then the mapping A i— A* defines an involution of the cone of dominant weights A"*" 
(2.7). The Macdonald polynomials labeled by A and A* are (for a,g real) related by 
complex conjugation 



(4.17) pa(x)=pv(x), AgA+. 

Indeed, the weight —A lies in the S'Ar+i-orbit of A* and the vector —p lies in the 
S'iv+i-orbit of p, from which it is concluded that 



mA(x) = mA(-x) = mA.(x), E+{p + A) = E;{p + A) = E+{p + A*). 

Combining this with the observation that /i* ^ X* if p ~< X then entails that the 
complex conjugate polynomial Pa(x) satisfies the same conditions of the type in 
(4.10a), (4.10b) as the Macdonald polynomial pa*(x), whence the equality in (4.17) 
follows by the uniqueness of the Macdonald polynomials. 

The upshot is that by passing to linear combinations of the form 

(4.18a) = (^^^)' ^^^M, 



(4.18b) *f = -V^ , AeA+ 



2i 

(thus selecting the real and imaginary parts of the wave functions \I'a (4.13)), we 
arrive at real-valued eigenfunctions for the discrete difference operators Hi, . . . , TYjv 
(3.2a), (3.9) 

(4.19) TY,^'^'/^ = Er{p + A)^'$'/^ A e A+ 

(r = 1, . . . ,N). Notice, however, that — in contrast to the complex wave functions 
(4.13) — the functions in (4.18a), (4.18b) do not diagonalize the operators Hf, . . . , 
(3.9) individually. When A runs through the alcove A^ (3.6), we of course count each 
eigenfunction ^£ and twice because 

(4.20) ^f* = *^, *f* = -*f, AeA+. 

Eliminating for this redundancy yields a real- valued orthonormal basis for the Hilbert 
space L^(p-|- A^i^), consisting of the wave functions with A e A^ mod * (i.e., 

we pick the weights from a fundamental domain in A^ with respect to the action of 
the involution *). The complex eigenbasis is recovered from the real eigenbasis by 
forming the combinations 

(4.21) = + AeA+. 
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5. Miscellanea 

5.1. The eigenfunction transform. The weight alcove K]^ (3.6) labels the eigen- 
basis (4.13). If we identify this weight alcove with the lattice p + then we are 
led to a Discrete Fourier- type Transformation in the Hilbert space LP'{p -\- Km), the 
kernel of which is determined by the cigenf unctions. 

Let : L^(p + A.'l^) — > L'^{p + ^Xi) discrete integral transformation 

(5.1a) (^/)(p+A)= HP + \P + P^)f{P + P^) 

with a kernel of the form 

(5.1b) T{p + X,p + p):^ *a(p + /^), 

where ^a(p+A*) is taken from (4.13). Furthermore, let £r '■ L^iP+^ti) ^ -^^(p+-Am) 
denote the multiplication operator 

(5.2) {£rf){p + p)^ Er{p + p)f{p + p) (r = 1, . . . , N) 

with Er{-) representing the real elementary symmetric function of (4.8). 

The main results of this paper may be conveniently summarized in the following 
three properties of the discrete integral transformation (5.1a), (5.1b) 

(5.3a) *:r = JT, JT* = JT-i 

and 

(5.3b) Tnr=ErT, r = l,...,Ar, 

where it is assumed that the parameters satisfy condition (3.5). The first property 
states that the transpose *JF of T is equal to JF, or in other words, that (the kernel 
of) the operator T is symmetric. This is a consequence of the symmetry relation 
in Proposition 4.6. The second property states that the adjoint T* of T in L'^{p-\- 
A^) equals the inverse of JF, or in other words, that the operator JF is unitary. 
This is a consequence of the orthonormality relations for the kernel ^a(p + p) in 
Proposition 4.8. Finally, the third property states that T simultaneously diagonalizes 
the discrete difference operators Hi, . . . ^Hn (3.1), (3.9) in L^(p-|- A^). This is seen 
by checking that both sides of (5.3b) act the same on the orthonormal eigenbasis ^a, 
A G A^ (cf. Proposition 4.7 and also Remark ii. at the end of Sect. 4). 

The map f ^ f '■— j determines a Discrete Fourier-type Transformation in 
L2(p + Aj^) of the form 

(5.4a) /(p-|-A) = (/,*a), AeA+ 

with the inversion formula given by 

(5.4b) /(p + //) = (/,*^), /.eA+. 

The coefficients j\ :— f{p + X), A G A^ solve the hnear interpolation problem of 
decomposing an arbitrary lattice function / e L^{p + A^) in terms of the wave 
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functions '^x, X E At, 



M 



(5.5) f^Yl A^A' /a = (/,*a). 



By passing to the real eigenbasis from Remark ii. at the end of Sect 4, which is 
given by (4.18a), (4.18b) with A e mod *, we arrive at analogs of the 

(discrete) Fourier cosine transform T'^ 

(5.6a) rip + A) = (/, ^r^), A e A+ mod * 

(5.6b) /(p + /i) = (r,^rj), /ieA+ mod* 

and Fourier sine transform 



(5.7a) /^(p + A) = (/,M/f), AgA+ mod 



* 



(5.7b) /(p + ^) = (/^M/^), ^eA+ mod*, 

respectively. Here it is assumed that / G L'^{p + A~lj) is a lattice function that is 
symmetric (for the cosine transform) or antisymmetric (for the sine transform) with 
respect to the action of the involution * on the lattice p + A^, respectively. (The 
involution * acts on the lattice p + Aj^ by (p + A) i— > (p + A*), cf. Remark ii. at the 
end of Sect. 4.) 

The integral transformation (5.1a), (5.1b) conjugates (cf. (5.3b)) the quan- 
tum Hamiltonians Hi, . . . , Hn (3.1), (3.9) to the multiplication operators £i, . . . ,£n 
(5.2). The spectrum of Hr in L^(p + Aj^) is (thus) given by the range of the real 
elementary symmetric function Er{-) — '^i.^Sn+iM ^^^^ lattice p + A~j^. 

The eigenfunction transform amounts to the quantum counterpart of the action- 
angle transformation 4> for the classical system, which was found in explicit form by 
Ruijsenaars [R3]. 

The action-angle transform at issue constitutes an (anti)symplectomorphism z z 
of the classical phase space (CP^, cvr), giving rise to new canonical coordinates (x, p) 
of the form (cf. (2.19a), (2.19b)) 

(5.8a) (a,-,p) = (27r/a - (iV + 1)^) , /^'^ ^ ,^ 



U^o H \- \zn\ 



(5.8b) e 





l^ol 







(j = 1, . . . , A^) on an open dense patch {[^o : ■ ■ ■ : z^] | 7^ 0, j = 0, . . . , A^} of 
CP^. The classical Hamiltonians ?ir(x, p) (2.15) become in these new coordinates 
of the form [R3] 

(5.9) ?i:^(x,p)= cosq;(i/,p), r = l,...,N. 

(This is to be compared with the jF-transformation of the quantum Hamiltonian to 
the multiplication operator 8^ in (5.3b).) That is, in the new coordinates the classical 
Hamiltonians depend only on the action variables p e and are independent of the 
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angle variables x e T = E/{2nQ) (cf. Sect. 2.3). Thus, the "spectrum" (i.e. the 
range) of the classical Hamiltonian Tir (2.15) on the compactified phase space CP'^ 
is given by the range of Ti^ (5.9) on the convex polytope {p | p € S^} (cf. (2.22)). 
(Notice that we may again continue the action variables p (5.8a) smoothly to the 
whole of CP^ unhke the angle variables x (5.8b), cf. Sect. 2.3.) At this point it 
is worthwhile to mention that the convexity of the range of the action variables 
for our model is in agreement with the general convexity results for Hamiltonian 
systems on compact symplectic manifolds due to Atiyah [At] and Guillemin-Sternberg 
[GSl]. We see that — roughly speaking — the quantization of the model discretizes the 
spectrum of the Hamiltonians as if the action variables p get localized on the lattice 
p + A^ <zT,g. It was furthermore shown by Ruijsenaars [R3], that the vertices p = p, 
p + Mur (r = 1, . . . , A^) of the convex polytope {p | p G S^} correspond to the 
equilibrium points of the flows generated by the classical Hamiltonians Tir (2.15) (cf. 
also Sect. 5.2 below). This state of affairs is again in agreement with the general 
picture presented by Atiyah and Guillemin-Sternberg [At, GSl]. 

A remarkable property of the action-angle transform is that it in fact defines an 
(anti) symplectic involution on {CF^,ur) (i.e. the model is "self-dual" in the termi- 
nology of Ruijsenaars) [R3]. The analogous property of the eigenfunction transform 
^ (5.1a), (5.1b) states that the corresponding discrete integral transform defines a 
(discrete) Fourier-type involution (i.e. an involution up to complex conjugation) in 
the Hilbert space L^(p -|- A^) 

(5.10) *:r*:r = id, 

cf. (5.3a). (The kernel of and are the same up to complex conjugation.) 
5.2. Ground-state vs nicLximal energy wave function. Let 

N+l 

(5.11) E{:x.) ^^cos{aXj). 

Notice that on the hyperplane Xi + - ■ --l-Xjv+i = this function coincides with the first 
(real) elementary symmetric function £'i(x) = cosq;(z/, x) (cf. (4.8)). It 

is not very difficult to infer that the critical points of E{'k) on the simplex (2.22) 
are located at the A^ -|- 1 vertices 

(5.12) p and p + Mur, r = l,...,N, 

where M = ^ — {N + l)g. The (critical) values of E{x) evaluated at these vertices 
read (cf. (4.6)) 



(5.13a) E{p) 



sinf (AT + l) 



sin(f) 
27rr 

We thus conclude (cf. Proposition 4.7) that the maximal eigenvalue of the Hamil- 
tonian Hi (3.1), (3.9) in the Hilbert space L^(p -|- A^) has the value -E'i(p) = 



(5.13b) Eip + Mur) = cos(^^)£;(p), r^l,...,N. 
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^™ ?' ^/Ien^^^ • The corresponding eigenfunction is given by the factorized wave function 

^0 in (4-4). For N odd the minimal eigenvalue reads Ei(p + Ma;(Ar+i)/2) = —Ei[p) 
whereas for even the minimal eigenvalue is twofold degenerate and given by 
E^{p + Mujm/2) = Ei{p + = cos(^)Ei(p). 

The "critical" or "vertex" eigenvalues of the Hamiltonian Tii in (5.13a) and (5.13b) 
coincide with the equilibrium values of the corresponding classical Hamiltonian 7ii 
(2.15) at the stationary points computed by Ruijsenaars [R3, Sect. 5.3] (cf. Sect. 5.1 
above). In particular, the global minimal and maximal energies of the Hamiltonian 
Til /Til read the same at the quantum level as they do at the classical level. In other 
words, the energy levels get discrctized at the quantum level (cf. Sect. 5.1 above) but 
there is no shift of the energy spectrum due to the quantization (as e.g. in the case 
of a harmonic oscillator). 

Prom a physical point of view it is often somewhat more natural to work with a 
nonnegative Hamiltonian. This can be achieved by passing to difference operators of 
the form 

(5.14) Hr = Er{p)-1ir 

2 ^ (k(x) + K(-x) 



- E 



r — 1, . . . ,N (where we have again employed the Macdonald identity from the proof 
of Proposition 4.3 to pass from the first to the second formula on the r.h.s.). The 
factorized eigenfunction ^'o (4.4) amounts to the ground-state wave function for the 
Hamiltonian Tii (5.14), with the corresponding eigenvalue being equal to zero. 



5.3. The two-particle solution. In the case of two particles, i.e. for N — 1, the 
quantum version of the compactified Ruijsenaars-Schneider model was introduced and 

solved already several years ago by Ruijsenaars in the survey paper [R2] (see Sect. 
3C2). It is instructive to view how, in this special situation, our results reproduce 
those previously obtained by Ruijsenaars. The difference operator Hi (3.1) (= 7if 
(3.2a) = (3.2b)) reduces to 

(5.15) H = ritt'T ^^ r°'tyY + 



sm[—) J \ sm(^ 2 



sin|(a; -5()y/2 _^ ^sin f (x + ^f) \ 1/2 

in(f) ) 



6 

sml'^) ) \ smi 



2 / """V 2 



with X — x\ — X2- Por 

277- 

(5.16) Q;,y>0 and 2^ = MeN\{0}, 

a. 
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the operator (5.15) becomes self-adjoint upon restriction to the Hilbert space 
LF'{g + A^) over the lattice 

(5.17) ^ + A+ ={^,^+1,^ + 2,... ,y + M} 

(cf. Sect. 3.2). An orthonormal basis of eigenfunctions '^i : [g + A^) — > R for 
(5.15) is given by (cf. Sect. 4, in particular Eqs. (4.13), (4.1) and Proposition 4.4) 

(5.18) ^>i{g + m) = _L^aV2(^) aV2(^) p^(cos '^{g + m)), 
Z,m = 0, . . . ,M, with 

A(m) = 



sinf(m + ff) (2g : sin^), 
sin(^) (1 : sina)„ 

sin f (m + 5') -Q sin ^{2g + j - 1) 



2 

2V"^"' ^' 1 I 2 ' 

sin(f) n sin(f) 



PKcos |(x)) = g^(--«)/2302 ( ^'''^f/ ^'""^ ; , g = e^" 



and 



M-l 



jVo = 2^(1 + g : sin,,)M-i = 2^^ J] o + 



fc=i 

Specifically, one has that (cf. Proposition 4.7) 

a 

(5.19) 7^^; = 2cos-(5 + 0^/, / = 0,...,M 
and that (cf. Proposition 4.8) 

M 

(5.20) J]^i(5 + m)*fe(^ + m) = (5,,fc, /, A; G {0, . . . , M}. 

m=0 

Clearly, the maximal and minimal eigenvalues of the Hamiltonian H. are given by 
2cos(^) and 2 cos f(^ + M) = -2cos(^), respectively (cf. Sect. 5.2). 

The (Ai-type Macdonald) polynomials Pi(cos ^{x)) coincide up to a normalization 
factor with the g-ultrasphcrical polynomials and can thus be explicitely written (in 
various ways) in terms of terminating basic hypergeometric series [GR, KS]. For our 
purposes it is convenient to employ the above representation in terms of a terminating 
302 series as this manifestly demonstrates the symmetry of the wave function '^i{g + 
m) (5.18) with respect to an interchange of / and m (cf. Proposition 4.6). The 
orthogonality relations (5.20) for the basis ^'i, / = 0, . . . , M boil down to reductions 
to special parameter values of well-known discrete orthogonality relations for the 
g-Racah polynomials due to Askey and Wilson [AS, GR, KS]. 

Notice that the rank-one case iV = 1 is very special in the sense that the reflection 
A I— > —A lies in the Weyl group (it amounts to the Weyl-permutation {xi,X2) ^ 
{x2.,xi) restricted to the hyperplane E — {{xi,X2) G | xi + 2:2 = 0}). As a 
consequence, the involution * of Remark ii at the end of Sect. 4 reduces in this 
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special case to the identity. Indeed, we have for = 1 that the wave function "ifx 
(4.13) is real-valued, whence = (4.18a) and = (4.18b). In other words, 
the "Discrete Fourier transform" = JF* (5.4a), (5.4b) and the "Discrete Fourier 
cosine transform" J-''^ (5.6a), (5.6b) coincide in the rank-one case and the "Discrete 
Fourier sine transform" (5.7a), (5.7b) collapses. 

5.4. Geometric quantization. In the light of the fact that the phase space for the 
classical compactified Ruijsenaars- Schneider model is given by the complex projective 
space CP^ equipped with the renormalized Fubini-Study symplectic form cur (2.20), 
it is natural to ask oneself the question as to what extent our results may be recovered 
within the realms of geometric quantization. In this formalism a Hilbert space is 
associated to the classical phase space {CF^,ujr) in two steps (see e.g. [Si, Hu]). In 
the first step (prequantization) , the question is to construct a Hermitian line bundle 
L with connection V over CP^ such that the curvature of V equals cur. Such a line 
bundle exists provided lur satisfies the integrality condition 



where the integration is over a complex projective line in CP^ (or more generally an 
integral two-cycle in the homology basis). (Geometrically, this condition means that 
uiR belongs to an integer cohomology class: [ujr] G if^(CP^,Z).) The prcquantum 
Hilbert space now consists of the space of sections of the line bundle L (where the 
measure of integration is taken to be the Liouville volume form associated to 0;^). 
After recaUing that the normalization of ur is such that § lor ~ At^B? (2.21), it is 
seen that the integrality condition in (5.21) amounts precisely to our quantization 
condition (3.5). (This observation was already made by Ruijsenaars in [R3, Sect. 



Unfortunately, the (prequantum) Hilbert space thus obtained is too big. Rough- 
ly speaking, it corresponds to an "L^ space over the phase space" whereas from a 
physical point of view one is interested rather in the analog of an "L^ space over 
the configuration space" . In the second step of the quantization procedure the pre- 
quantum Hilbert space has to be downsized so as to produce the physical Hilbert 
space. To this end it is needed to exploit the fact that ,ujr) is a Kahler mani- 
fold and, as such, carries a natural Kahler polarization. Specifically, as the physical 
Hilbert space one picks the subspace of the prequantum Hilbert space consisting of 
all sections of the line bundle L that arc covariantly constant (with respect to 
the connection V) along the leaves of the (standard) Kahler polarization on CP^. 
The result is a physical Hilbert space li-hoi consisting of the holomorphic sections of 
L. The dimension of the space of holomorphic sections Tihoi follows from a classical 
result (viz. a Riemann-Roch-type formula) due to Hirzebruch and Kodaira [HK] 



(5.21) 




1.3].) 



(5.22) 
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which corresponds nicely to the dimension of our Hilbert space L^(p + Aj^) in (3.8). 
It is possible to realize the Hilbert space Ti-hoi more explicitly, as the space of holomor- 
phic sections may be identified with the space of functions of the form [HK, Si, Hu] 

^ ' {i + \zi\'+---+\zN?r' 

where p{zi, . . . ,Zn) denotes an arbitrary polynomial of degree at most M in the 
affine CP^ coordinates {zi, - ■ ■ , Z]si). In this representation the integration of the 
inner product is with respect to the volume form 

(5.24) (1 + H h I^Ar^-^^+^^rf^i ^dzl^■■■ ^dzN ^ dzN- 

It would be very interesting to extend the analysis further so as to include a de- 
scription of the quantum Hamiltonians and their eigenfunctions within the framework 
of geometric quantization and to compare the results with the approach taken in the 
present paper. In this connection it is expected that our lattice p + A^^ may be re- 
covered geometrically as the so-called Bohr-Sommerfeld set — see Guillcmin-Sternberg 
[GS2] — associated to the symplectic embedding of x T into CP^ given by (2.18). 
Furthermore, the embedding in question induces a "real polarization with singulari- 
ties" on CP^ (cf. [GS2]). (This real polarization becomes singular at the boundary 
hyperplanes Zj = of the open dense patch {[zq : ■ ■ ■ : z^] | 2;^ 7^ 0, j = 0, . . . , A^} C 
CP^.) The above-mentioned correspondence between (the dimensions of) the Hilbert 
space L^(p + K\j) and the Hilbert space of holomorphic sections of the line bundle L 
(viz. l-Lhoi) suggests that (the dimension of) the Hilbert space associated to (CP^, ujr) 
via geometric quantization does not depend on the choice of polarization to be either 
the standard Kahler polarization or the "real polarization with singularities" stem- 
ming from the embedding x T CP^. This is in correspondence with the more 
general "invariance of polarization" results for the geometric quantization of complex 
flag manifolds due to Guillemin and Sternberg [GS2] . 

5.5. Connections to integrable field theories. The compactified Ruijsenaars- 
Schneider model is related in various ways to well-known infinite-dimensional inte- 
grable systems. For instance, in [R5] it was shown that at the classical level the 
(r-functions of) single-solitons for the Kadomtscv-Pctviashvilli and 2D Toda hierar- 
chy (cf. [DKJM, JM, Hoi, Ho2, ZC]) may be described in terms of the equilibrium 
behavior of the classical compactified Ruijsenaars-Schneider molecule with an appro- 
priate center-of-mass motion. Multi-solitons arise in this picture — in a nutshell — by 
passing to composite integrable Ruijsenaars-Schneider-type (mi + ■ ■ ■ + m„)-particle 
systems that arc built of n (the number of solitons) interacting compactified trigono- 
metric Ruijsenaars-Schneider molecules in their ground state [R5]. 

In [GN] it was furthermore argued that formally the (quantum) compactified 
Ruijsenaars-Schneider model can be obtained by Hamiltonian reduction from an 
infinite-dimensional system on the cotangent bundle over a central extension of the 
loop group SU (A" -|- 1). The latter paper also indicates some intriguing relations be- 
tween, on the one hand, the quantum compactified Ruijsenaars-Schneider model and, 
on the other hand, a gauged SU{N + 1)/ SU{N +1) Wess-Zumino-Witten topological 



30 



J. F. VAN DIEJEN AND L. VINET 



quantum field theory on a cylinder and a Chern- Simons theory with gauge group 
SU{N + 1) on a three-fold that is the product of an interval and a real two-torus. 

5.6. Bispectrality. The symmetry of the wave function '^x{p + lA (4-13) with re- 
spect to an interchange of A and ix (cf. Proposition 4.6) has as consequence that 
it satisfies the same discrete difference eigenvalue equations in the "spectral" vari- 
able A as it does in the "spatial" variable ii (cf. Proposition 4.7 and the proof of 
Proposition 4.8). In other words, we are dealing with a multivariate doubly discrete 
finite-dimensional bispectral problem in the sense of Duistermaat and Griinbaum 
[DG, W, G]. 

Appendix A. Truncated and terminating Aomoto-Ito-Macdonald sums 

In this appendix a finitely truncated version of a recent summation formula due to 
Aomoto, Ito, and Macdonald [Ao, I, M5] (see also [Ka]) is derived. In Sect. 4 we used 
this finite Aomoto-Ito-Macdonald-type sum to arrive at a compact product formula 
for the normalization constant of our factorizcd wave function \l/o (4.4) (cf. Proposi- 
tion 4.4). When formulating the summation formulas in question it is convenient to 
employ the g-shifted factorial defined by (see e.g. [GR]) 

r (1 - aq){l - aq^) ■■■{!- aq'^''^), m = 1, 2, 3, . . . 
(a;?)m= < 1> ^ "^ = 

I (l-oq-l)(l-aq-2)-(l-oq'") ' m = -1, -2, -3, . . . 

(where for negative m it is assumed that a,q & C are such that the denominators do 
not vanish), 

oo 

(a;g)oo = n(l-0' (0< |g| <1) 

n=l 

and 

(ai, . . . , Uk; q)m = (oi; q)m--- {ak; q)m, me ZD {oo}. 

Proposition A.l (The Aomoto-Ito-Macdonald sum [Ao, I, M5]). For < q < 1, 
Re{g) <0 andze C^+^ with g+ {a,z) + for all a e A+, one has that 

(where the sum is taken over all weights in A (2.5) and the vector p is given by 
(2.13);, with 

e(z) = n ^S^v = C, gC-^ g)oo 



aeA+ 



and 



7 



^11 (gl-(a,p)^^-(a,p>.^^^ 



a€A+ 
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where Sg, :—lifa is simple (of. (2.2) ) and 5^ '■— otherwise. Furthermore, the series 
on the l.h.s. converges in absolute value. 

The conditions on q, g ensure that the series on the l.h.s. converges absolutely and 

the genericity restrictions on z guarantee that all denominators are nonzero. The 
sum of Proposition A.l was first considered by Aomoto [Ao], who showed that it 
can be evaluated as a product of the quasi-periodic factor 0(z) and a z-independent 
constant. The value of this constant (viz. 7), was subsequently conjectured by Ito [I]. 
In its present form the statement of the above proposition is due to Macdonald [M5], 
who derived it by linking a constant term identity of Cherednik [CI] to a generalized 
Poincare series type formula for affine Weyl groups due to Matsumoto [Ma] . 

After division of both sides of the Aomoto-Ito- Macdonald summation formula by 
the middle term on the l.h.s. corresponding to /i — 0, one arrives at the identity 



(a,/i> 



7 



n 



aeA+ 



(^l+(a,z)^^l-(a,z).g)^ 
(gl-c/+(a,z)^ gl-g-(a,z). 



where it is assumed that < q < 1, Re{g) < and that z G C^"*"^ satisfies the 
genericity conditions (a, z) ^ jj^^Z, and g + (a.z) — 1 <^ N + j^^Z, for aU a G Ajq 
(to ensure that there is no division by zero). It is instructive to observe that the 
proportionality constant 7 (see above) on the r.h.s. may be rewritten in a somewhat 
more compact (but less elegant) form by canceling common factors in the numerator 
and denominator: 



(A-2) ^=(^+1)7;^^ n 



We will now show that by specializing the vector z to the value p (2.13), the sum 
in (A.l) over the weight lattice A (2.5) truncates to a sum over the dominant cone 
A+ (2.7). 

Proposition A. 2 (A Truncated Aomoto-Ito-Macdonald sum). Let < g < 1 and 
Re{g) < such that g — (a, p) is not a positive integer modulo for all a G A^. 

Then 



n/ l+(a,p) S^+g-(a,p) \ 
(^-(a,p>^^l-S+(a,p>. 



aeA+ 



>;?)c 

^ (-^^) n ^ 

and the series on the l.h.s. converges in absolute value. 
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Proof. The conditions on g (and q) ensure that after substituting z — p (2.13) in (A.l) 

all terms remain finite and the series converges in absolute value (as a consequence of 
Proposition A.l). The resulting series on the l.h.s. now truncates because the terms 
become zero for e A \ A"*". This is because for e A \ A"*" there exists a simple root 
3Lj (2.2) for which (a^, /x) is a negative integer and hence we pick up a zero from the 
factor 



(which is zero for (a^-, /x) < 0). The expressions for the r.h.s. are obtained from that of 
(A.l) (with z = p) by canceling common factors in numerator and denominator. □ 

A further reduction arises when we specialize the parameters q and g in such a way 
that = 1 for some positive integer M. The sum of Proposition A. 2 then 

terminates to a sum over the integral alcove A^ (3.6). 

Proposition A. 3 (A terminating Aomoto-Ito-Macdonald sum). Let 

2m 

' ^ ^"P^ ff(Ar + i) + M ^ 
with g > and M a positive integer. Then 

^(a,p+M>^ (g^+<^'^^g)(a,^) 



/ 1 _ n\^'P+M \ 



= (iv + i) n s)"-!- 

l<n<N 

Proof. Let us first substitute 
(A-3) 9 = -^r-- + 



N + l {N+l) log(g) 

with < g < 1 and M a positive integer, in the summation formula of Proposi- 
tion A. 2. Notice that this value of g satisfies both the convergence criterion Re{g) < 
as well as the regularity condition that — (a, p) — 1 ^ N + tj^^Z for all a e A^. 

The sum over the dominant cone A"*" (2.7) then terminates to a sum over the integral 
alcove A^ (3.6) because all terms become zero for ji G A+ \ A^j. Indeed, for the 
above value of g we have that = 1, so we pick up a zero from the factor 

W lH){a.max,IJ.) W ) yy'Ml-/*JV+l 

when (a^ax, A*) = A*i — A*Ar+i > M. To arrive at the expression for the r.h.s. one uses 
that 

(g-"^;g)oo " (g(^+i-")^+^g)oo " 

l<n<Ar i^Joo i<n<N J'i/oo i<„<7v 

We thus see that the terminating sum of the proposition holds for q given by 
ew{ ^N+i)l+M ) wi^^ ^e(f) = -M/{N + 1) and Im{g) < (solve (A.3) for q). By 
exploiting the analyticity in g it is possible to extend the terminating sum to generic 
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complex g. The restriction to positive real values of g ensures that all numerators 
and denominators are nonzero. □ 

In trigonometric notation with q — e*", the summation formula of Proposition A. 3 
becomes 



(A-*) E n 



2 P + n 2 ^ ^ ^ ~ 

m=l 

a6A+ sin(|(a, p)) H sin |((a, p) - <7 + m) 

m=l 

= 2^(^-^) (N + l) II sin I (m + ng) 

l<m<M-l 

l<n<N 

for a = 27t/ {g{N + 1) + M) with M a positive integer and g > 0. (Here we have used 
the convention that empty products are equal to 1.) This is precisely the summation 
formula of Proposition 4.4. 

Remarks: i. For = 1 the Aomoto-Ito-Macdonald-type sums of Eq. (A.l), Propo- 
sition A. 2 and Proposition A. 3 reduce to 

(with < g < 1, Re{g) < and z ^ -^^Z, g±z-l^N+ j^^Z), 

.A 6) - (g^^; q)m _ ^ {q'^<^; g)oo 

(with < g < 1, Re{g) < 0) and 



m=0 



(with g = 9 respectively. The sums in (A.5), (A.6) and (A.7) 

are well-poised 2'02, 2<f>i and terminating 20i sums that arise as reductions of Bailey's 
QipQ, Rogers' 605 and Rogers' terminating very- well-poised sums, respectively 
[GR]. 

ii. The summation formula of Aomoto, Ito and Macdonald in [Ao, I, M5] is in 
fact more general than as stated in Proposition A.l. This is because they consider 
sums associated to an arbitrary reduced integral root system. The formulation of 
Proposition A.l corresponds to the restriction of [Ao, I, M5] to case of the Aj^f series. 
In [D3] a further extension of the Aomoto-Ito- Macdonald sum for the nonreduced BC- 
type root systems was studied together with corresponding truncated and terminating 
variants. 
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Appendix B. Bilinear summation identities for Macdonald's 

symmetric functions 

The purpose of this appendix is twofold. Firstly, it serves to collect some basic 
facts on the Macdonald symmetric functions that were needed in Sect 4. For a more 
complete treatment of this material and proofs the reader is referred to [M4, Ch. 
6]. Secondly, the appendix allows us to reformulate some of our results from the 
perspective of algebraic combinatorics. This gives rise to a novel system of bilinear 
summation identities for the Macdonald symmetric functions (cf. Proposition B.2 
below) . 

Let 

(B.l) Dr^t^^-^y^ Yl r = l,...,7V + l, 

Jc{i,...,Af+i} jeJ ^ ' 
\J\=r 

where \J\ denotes the cardinality of the index set J C {1,... ,A^ + 1} and Tjq = 

{Tj,qf){Zl, . . . , Zn+i) = f{Zi, . . . , Zj-i, qZj, Zj+i, . . . , zn+i)- 

Let n = (ni, n2, . . . , un+i) G N^"*"^ be a partition, i.e., let the components (or parts) 
be ordered as ni > 712 > • • • > un+i > 0. The monomial symmetric function mn(z) 
associated to n is then defined as 

(B.2) m„(z)= Yl ^r---4:r, 

me5Ar+i(n) 

where the sum is over the orbit of n under the action of the permutation group Sn+i 
on the components. The basis of monomial symmetric functions {nrin} inherits a 
partial order from the dominance partial order of the partitions defined by 

(B.3) m n iff |m| = |n| and mi H h < rii H \- rik 

ior k — 1, . . . ,N, where |n| :— rii + • ■ ■ + n^v+i denotes the weight of the partition. 

Proposition B.l (Triangularity [M4]). The q-difference operators Di, ... , L'jv+i 
commute and are triangular with respect to the basis of monomial symmetric functions 

DrTUn = Y [-Drln.mrUm, [-Dr]n,m £ Q{q, t). 

m^n 

Furthermore, the diagonal matrix elements (eigenvalues) [-D,]n,n o'^e given by 

Jc{i,...,N+i} jeJ 

For a quick proof of the polynomiality of (ZJ^^TT-n) (z) in z one may use that this 
rational expression is regular as a function of z e C^"*"^ due to the permutation 
symmetry. The triangular form of the monomial expansion for (iD^mn) (z) and the 
diagonal matrix elements [Dr]n,n then follow from the asymptotics for z to infinity. 



THE COMPACT QUANTUM RUIJSENAARS-SCHNEIDER MODEL 35 

In the simplest situation, i.e. for n = (so rrin — 1), the monomial expansion 
of Proposition B.l reduces to the (AN-type) Macdonald identity (cf. [Ml, Theorem 
(2.8)1) 

(B.4) y: n(|ff)= E n 

Jc{i-... ,N+i} jeJ ^ ^ ^ Jc{i,... ,iv+i} jeJ 

,/|=r fe^J \j\=r 

with r = 1, . . . , + 1. Substitution of t = e*"^ and Zj = e'"^\ j = 1,... ,N + 1 

rN 

produces upon division by t~ the Macdonald identity used in the proof of Propo- 
sition 4.3 for r = 1, . . . ,N (the {N + l)-th identity in (B.4) is in fact trivial, as in 
that case the product on the l.h.s. becomes empty and one ends up with the equality 

^N{N+l)/2 ^ fN^N-l . . .f . l^j 

The Macdonald symmetric functions arc now defined as the (joint) eigenfunctions 
of the commuting operators Di, . . . , D^^i. Such a definition makes sense because 
the eigenvalues En,r{Q,'t) nondegenerate (and hence semisimple) over the field 



(B.5a) p„(r) = J] 



Definition (Macdonald symmetric functions). For a partition n in N-^^^, the Mac- 
donald symmetric function Pn(z) is the symmetric polynomial of the form 

(i) p„(z) = m„(z) + ^ Cn,m{q,t) m^{z) with c^^^{q,t) e Q{q,t) 

m-<n 

such that 

(n) DrP^ = E^Aq, t)p^, r^l,... ,N + 1. 

Two important properties of the Macdonald symmetric functions are the evaluation 
formula (also referred to as the specialization formula) [M4, Ch. 6: Eqs. (6.11), 
(6.11')] 

11 (fk-j-n) 

l<j<k<N+l ^^)rij-nk 

with r = {t'^ ,t^~^, ... ,t,l) and the symmetry relation [M4, Ch. 6: Eq. (6.6)] 

(B.5b) Pm(rg") = P„(Tg™) 

with 

(B.6) Pn(z) = P„(z)/p„(t) 

and rg" = {t^q''\t^-'^q''\ . . . , ig"^, g"^+i). 

The Macdonald symmetric function p„(z) is homogeneous of degree |n| in z and 

P„+ei+...+e^+i(z) = (^1 • • • Zn+i)pM- 

Projection to a homogeneous function of degree zero in z 
(B.7a) Pa(z) = (^i---z^+0-l"l/(^+^)p„(z), 

(B.7b) A = n--J^(ei + --- + ejv+i) 
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gives rise to the Macdonald polynomials p^, A e A+ (2.7) associated to the root system 
A]\[ [M2, M3]. The functions Pa(z) (B.7a), (B.7b) are related to the trigonometric 
Macdonald polynomials px{x.) of Sect. 4.2 via the trigonometric substitution 

(B.8a) t = q = e^" 

(B.8b) Zj = e^"^^' , J = 1, . . . , iV + 1. 

More specifically, by substituting (B.8a), (B.Sb) the functions Pa(z) (B.7a), (B.7b) 
pass over to trigonometric polynomials Pa(x) of the form in (4.10a), (4.10b). (The 
difference equations in (4.10b) are equivalent to the g-difference equations for Pa(z), 
originating from the g-difference equations D^p^ — En,riQ,t) in the above def- 
inition of the Macdonald symmetric function Pn(z), upon substitution of (B.8a), 
(B.Sb)). The symmetry relations in (4.12) for the renormalized trigonometric Mac- 
donald polynomials -Pa(x) (4-11) are an immediate consequence of evaluation formula 
(B.5a) and the symmetry relation (B.Sb). The real-valuedness of the expansion coef- 
ficients CA,ji( in (4.10a) follows from the fact that the Macdonald symmetric functions 
Pn(z) are invariant with respect to the parameter inversion {q,t) {q~^,t~^), see 
[M4, Ch 6: Eq. (4.14) (iv)]. 

Translating back the orthogonality relations from Sect. 4.2 leads us to the following 
system of bilinear summation identities for the Macdonald symmetric functions Pn(z). 

Proposition B.2 (Bihnear summation identities). Let n, k be partitions in 
with rii — riN+i, ki — kN+i < M e N \ {0}. Then we have that 

X > |m|(|n|-|k|) , 



2^ 

M>mj >■■■> 



q ^+^p„(rg'")pk(r-'g-'")A(m) 

i/ k 7^ n mod Z (ei + ■ ■ ■ + Gat+i) 



JV+1=" 

^ t " 2 ' "A/'oA/'(n) if k = n mod Z{ei + h Gat+i) 

as rational identity in q^^ and t subject to relation t^^^q^ — 1, where 

.±l„±m 



and 



j\f{n) = n 



l<i<ifc<Ar+l ,V)mj-mk 

{t^+^-i,qt^-i-\q)n^.n, 
{t^-^,qt^-i;q)n.^. 



l<j<k<N+l ^ ' ^ ' ^^"J' 

Afo = {N+1) H {qt^;q)M-i. 

(For k = n mod Z (ei + ■ ■ ■ + Gat+i) the relation is actually rational in q itself.) 

Proof. First the projection formulas in (B.7a), (B.7b) are used to rewrite the stated 
bilinear summation identities in terms of the A^v-type Macdonald polynomials Pa, 
A e A'Ij. The trigonometric substitution in (B.8a), (B.Sb) then leads us back to the 
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discrete orthogonality relations (4.15a) for the trigonometric Macdonald polynomials 

Pxi'x), A G of Sect. 4.2. This proves our bilinear summation formulas for the 
Macdonald symmetric functions Pn(z) for t = q^, q = Gxp{jj^^^^) with g > 0. 
Analytic continuation in g entails the formulation of the proposition. □ 

If we specialize the formula of Proposition B.2 to the case that n = m = 0, then we 
arrive at the following rational identity in q,t subject to the relation t^^^q'^ = 1 

(B.9) Yl ^("i) = -^0' 



M>mi >--->m 



N+l = 



which amounts to the terminating Aomoto-Ito-Macdonald sum of Proposition A. 3. 

The complex conjugation relation in (4.17), for the trigonometric Macdonald poly- 
nomials Pa(x) of Sect. 4.2, translates to a corresponding relation for the Macdonald 
symmetric functions Pn(z) defined above. Specifically, if we associate to a partition 
n e N^"*"^ a contragredient partition n* e with parts given by 

(B.IO) n* = ni - n(N+2~j), j = 1, . . . , + 1 

(see Figure 4), then we have that 

(B.ll) p„, (z) = {z,--- ^^+0^ p„(z-^) = n, 

where z~^ := (zi^,... ,^^^^). Observe that the mapping n i— n* is involutive 
modulo Z(ei + • • • + Gn+i), i.e., the contragredient partition of n* is equal to n up 

to a possible integer multiple of the vector ei + h stv+i (in case un+i > 0). The 

verification of (B.ll) goes along lines very similar to the proof of Proposition B.2. 
First one uses the projective relation between the Macdonald symmetric functions 
and the Ajy-type Macdonald polynomials Pa in (B.7a), (B.7a), to conclude that 
after the trigonometric substitution (B.8a), (B.8b) the relation in (B.ll) reduces to 
(4.17). Notice to this end that if A is the projection of a partition n e N^^'^ C ]R^+^ 
on the hyperplanc E (2.1) (cf. (B.7b)), then A* (4.16b) amounts to the projection of 
the contragredient partition n* (B.IO) onto E. This proves Eq. (B.ll) for t = q^ with 
g > and q, zj on the unit circle. Analytic continuation then entails that Eq. (B.ll) 
holds identically as an equality that is polynomial in z and rational in q, t. 

For Schur functions (t — q) the formula in (B.ll) is well-known, sec e.g. Stanley 
[St, Eq. (11)]; it expresses an equivalence between (the characters of) the irreducible 
representation of SL{N + 1, C) associated to the partition n* and the representation 
contragredient to the one associated to n. We have not been able to locate a reference 
for the property (B.ll) applying to the general (g, t)-Macdonald symmetric functions, 
but most likely it was known already for this case too. 

With the aid of (B.ll), one rewrites the equality of Proposition B.2 in the form 

(B.12a) Yl ?-^'^^P„(rg'")PkN")A(m) 



M>m]^>--->mjY -1-1=0 



if ky^n* mod Z(eiH he^v+i) 

t^^'"2"^'''"A/'o A/'(n) if k = n* mod Z (ei H he^v+i), 
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Figure 4. The contragredient partition n* of a partition n e N-^+^. 

or equivalently in the normahzation of (B.6) (cf. also (4.15b)) 
(B.12b) g- '""^1"' P„(T?™)Pk(T?™) A(m) 

M>mi >--->mjy^_|_j=0 

_ r if ky^n* mod Z(eiH h eiv+i) 

i/k = n* mod Z(ei + --- + ejv+i), 

as a rational identity in t and (or even g if k* = n mod Z (ei + ■ ■ • + eAr+i)) 

subject to the relation t^^^q^ — 1, where n and k are partitions in with ni — 

riN+i, ki — k^^i < M. (To derive (B.12b) from (B.12a) one divides hy p^{t)p^{t) and 
uses that t-^l"l/2p^(r) = r^l^'l/^p^. (r) and that t^l"lAA(n)/(p„(r))2 = 1/A(n).) 

Remarks: i. In full generality, Macdonald defined his symmetric polynomials 
for an arbitrary integral root system [M3]. For the nonreduced BC root systems, 
Koornwinder [Ko2] subsequently found a further generalization leading to a class of 
Askey- Wilson polynomials [GR, KS] in several variables that contains all Macdonald 
polynomials associated to the classical root systems as special cases (cf. [Dl, Sect. 
5]). In [DS] finite-dimensional discrete orthogonality properties of a type analogous 
to those described by Proposition B.2 were studied for Koornwinder's generahzed 
BC Askcy-Wilson-Macdonald polynomials. In the case of the discrete orthogonality 
structure, the BC polynomials in question may be viewed as a multivariate gener- 
alization of the well-known g-Racah polynomials introduced by Askey and Wilson 
[AS, GR, KS]. 

ii. For t = q^ with g a nonnegative integer, the condition t^^^q^ = 1 implies that 
g is a root of unity. In this special case the properties of the Macdonald polynomials 
were studied by Kirillov, Jr. and Gherednik [Ki, G2]. In particular, Kirillov, Jr. 
connects the Macdonald polynomials at issue with the representation theory of the 
quantum group (quantized enveloping algebra) Uq{slN+i) for q a root of unity. 
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